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ABSTRACT 

In accretion-based models for Sgr A* the X-ray, infrared, and miUimeter emis- 
sion arise in a hot, geometricaUy thick accretion flow close to the black hole. The 
spectrum and size of the source depend on the black hole mass accretion rate 
M. Since Gillessen et al. have recently discovered a cloud moving toward Sgr A* 
that will arrive in summer 2013, M may increase from its present value Mq. 
We therefore reconsider the "best-bet" accretion model of Moscibrodzka et al., 
which is based on a general relativistic MHD flow model and fully relativistic 
radiative transfer, for a range of M. We flnd that for modest increases in M the 
characteristic ring of emission due to the photon orbit becomes brighter, more 
extended, and easier to detect by the planned Event Horizon Telescope submm 
VLBI experiment. If M > 8Mq this "silhouette" of the black hole will be hid- 
den beneath the synchrotron photosphere at 230GHz, and for M > WMq the 
silhouette is hidden at 345GHz. We also flnd that for M > 2Mo the near-horizon 
accretion flow becomes a persistent X-ray and mid-infrared source, and in the 
near-infrared Sgr A* will acquire a persistent component that is brighter than 
currently observed flares. 



1. Introduction 



^The recent discovery of a cloud moving towards the Galactic Center ( IGillessen et al. 

20121) creates a potential opportunity for testing models of Sgr A*. The > 3M^ cloud will 
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interact strongly with gas near nominal pericenter at 300AU 8000GM/c^ (M = 

black hole mass), and may change the black hole accretion rate M. Since the structure of 
the cloud and the surrounding medium are uncertain, possible outcomes range from very 
small changes in the accretion rate over timescales of decades to rapid, large changes in the 
accretion rate. 

The dynamical timescale at is td = (r^/(GM))-^/^ = 0.5yr, and the viscous timescale 
tvis = {oi^)~^{R/ HY ^ lOOyr ^ tpho^ assuming a = 0.05 and H/R = 0.3, i.e. a hot, radia- 
tively inefficient accretion flow. After an initial transient phase while the flow circularizes — 
accompanied by transient emission — it is natural to think the flow will settle into a steady 
state. The settling timescale could be as little as a few td^ and so the steady state may arrive 
as soon as mid-2014. If the resulting flow can be modeled as a steady disk, the excess mass 
will drain away on the viscous timescale, i.e. the source will remain bright well into the 
22nd century. It is therefore interesting to ask how changes in M will manifest themselves 
observationally. 

Current ob servations of Sgr A* show Fjy = 0.5 — 1 Jy at 1-50 GHz with a nearly flat 
spectral slope (IFalcke et al.l Il998l ): Fj^ ^ ^p=o.i7-o.3 The spectral slope becomes flatter 
and variable at 230-690 GHz, with p = —0.46 — 0.08, iMarrond l2006l ) . which is commonly 
interpreted as signaling a transition from optically thick to optically thin synchrotron emis- 
sion. The discovery of polarized emission (polarization fraction at level of a few per cent) 
at A = 1.3mm and subsequent measurement of Faraday rotation imply a model dependent 
limit 2 X 10"^ > M > 2 x lO~^M0yr~^ in the inner accretion flow (IBower et al.l [20051, 



Marrone et al.ll2006l ). Sgr A* is resolved by VLBI at 1.3mm where interste llar scattering is 
comparable to the intrinsic source FWHM 37l}^//as flDoeleman et al.ll2008l ). This is smaller 
than the apparent diameter of the event horizon ^ 55uas. Sgr A* fluctu ates rapidly in 



the near-infrared (NIR), with average Fj^ ^ ImJy (iDodds-Eden et al.ll201lh . It is not yet 
detected at mid-infrared (MIR) wavelengths (A < 8.6/xm, e.g. ISchodel et al.ll201ll ). In the 
X-ray Sgr A* exhibits flares with characteristic duration of order an hour and a duty cycle 
of order 5%. The up per limit for quiescent state X-ray emission is uLj, < 2.4 x lO^^ergs"^ 
(IBaganoff et al.ll2003l ). We are not aware of any secular trends in these observed properties 
of Sgr A*. 

O ne model of Sgr A* that flt s most observational constraints is our relativistic accretion 
model (IMoscibrodzka et al.ll2009l ). where submillimeter, IR, and X-ray emission arise in an 
optically thin, geometrically thick disk close to the event horizon; rad io emission is ass umed 
to arise nonthermally in a synchrotron photosphere at larger radius (lOzel et al.ll2000li ). but 
is not predicted by our model, which focuses on the inner parts of the accretion flow. The 
underlying flow model is a general relativistic magnetohydrodynamic (GRMHD) simulation 
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(IGammie et al.l 120031 . iNoble et all 120061. IShio 



calculated using Monte Carlo (IDolence et al. 



kawa 



et al.l I2OIII ) . The emerging radiation is 



20091 ) and ray- tracing schemes. The 1.3mm 



flux originates as Doppler boosted synchrotron emission from the approaching side of the 
disk between ^ IQGM/c^ and the innermost stable circular orbit (ISCO). 

Our relativistic disk model assumes a thermal electron distribution function, neglects 
thermal conduction, and assumes a constant ratio of ion to electron temperature T^/Tg. The 
model does not produce the observed IR flaring at the correct amplitud e, but can w ith the 
addition of a small nonthermal component in the distribution f unction (Leung 201ol ). If we 
flx the black hole ma ss {Mbh = 4.5 x lO^M©, Ichez et al.lboosh and distance {D = 8.4 kpc, 
Gillessen et al.ll2008l ). the remaining model parameters are the source inclination i, black 
hole spin (0 < a* < 1), T^/Tg, and M. We flx M so that the 1.3mm flux matches the 
observed 3 Jy. 

The relativistic accretion models are not tightly constrained by the data, but they 
reveal the following: (1) face-on models that reproduce the millimeter flux would look like 
rings and therefore, in VLBI data, have dips in visibilities on flxed intermediate baselines, 
while existing ob servations suggest that the ring; radius would need to vary to flt the data 



( IFish et al.l 1201 ll : see also the discussion of iBroderick et al.l I2OIII ) . More nearly edge-on 
models are therefore favored; (2) models with a* ^ 0.98 that reproduce the millimeter flux 
have a hot, dense inner disk that would overproduce X-rays via inverse Compton scattering. 
Lower spin, a* ^ 0.9, models are therefore favored; (3) the observed source size and flux 
flx the temperature of the emitting electrons Tg = Fj^c^ / {Ankiy^a'^) {a = the RMS size of 
the s ource on the sky) and this turns out to favor T^/Tg 3 models. The "best-bet" model 
from lMoscibrodzka et all koO^ has a* 0.94, i = 85deg, M = Mq - 2 x lO"^M0yr"^ and 



T,/Tg = 3. 

The relativistic disk model uses self-consistent dynamics and radiative transfer but is 
not unique. The electron distribution function is particularly poorly constrained. It is 
likely anisotropic, may contain multiple temperature components ( IRiquelme et al.ll2012l ) and 
power-law components, and vary in basic functional form with time and po sition. Alterna- 
tive accretion models (e.g. IBroderick et al.l I2OIII . IShcherbakov et al.l I2OIOI ) make different 
assumptions about the flow and/or distribution function and favor slightly different M,a*, 
and i. These models may respond differently to an increase in M. 

Indeed, radically different models may also flt the data. The model dynamics depends 
on the initial magnetic fleld distribution, particularly the distribu tion of vertical flux thro ugh 
the accretion disk; models with large vertical magnetic flux (e.g. iMcKinney et al.l 120121 ) are 
l ikely to respond quite differently to variations in the mass flux. A lso, jet models for Sgr A* 



(IFalcke fc Marko 



2000 



Loeb fc Waxmanll2007l . iFalcke et al.ll2009l ) posit a luminous jet and 
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a comparatively dim accretion disk. Again, these may respond differently to an influx of 
mass. 

How, then, does our relativistic disk model respond to an increase in Ml In this Letter 
we calculate the 1.3 and 0.87mm flux and source size as well as the spectrum that would 
result for the best-bet model over a range of M . One key question we seek to answer is 
whether a small increase in M would hide the event horizon (and the signature ring-like 
appearance of the photon orbit, also known as the shadow or silhouette of the event horizon) 
underneath a synchrotron photosphere. This mi^ht prevent det ection of the photon orbit 



by the planned Event Horizon Telescope ( IDoeleman et al.ll2008l ). Another key question is 



whether the increased M would make Sgr A* detectable in its quiescent state in the IR and 
X-ray. Below we describe variation of the flux and source morphology at 230GHz (1.3mm) 
and 345GHz (0.87mm) (§ [2j), describe variation of the spectrum (§ [3j), and flnally discuss 
which features of the results are likely to be most robust (§ [4]). 



2. Change of Sgr A* sub-mm luminosity and size for enhanced M 

How do we naively expect the mm disk image size and flux to respond to changes in 
Ml In our model, 1.3mm emission in Sgr A* is thermal synchrotron emission from plasma 
with optical depth ^ 1, near the ISCO. The true electron distribution undoubtedly 
contains nonthermal components (se e, Riquelme et al. 2012). Models with thermal + 



power-law distribution functions (e.g. iBroderick et al.ll201ll ) contain an 0(1/3) nonthermal 



contribution to the flux at 1.3mm, which hints at how uncertainty in the distribution function 
translates into uncertainties in the spectrum. 

The thermal synchrotron absorptivity is a^^a = ju/By^ where jy = ^^^^(^"^^^ + 
211/12^1/6^2 g^p(_^i/3)^ X = jy/jy^^ = 2 / 9 {cB / 271771^0)91 siu 9 , ^ is an angle between the 
magn etic field vector an d emitted photon, K2 is a modified Bessel function of the second 



kind ( iLeung et al.ll201ll ) and 2z/^0eme. Near 1.3mm, X ^ 1, and the emissivity is 

nearly independent of frequency, so oc u^rieB and a^y^a oc iy~'^neBQ~^. 

We will assume that oc Mr~^/'^^ Gg oc 1/r, and /3 ^ const., so that B oc M'^/'^r~^/^ for 
r > GM/c^^ and ignore relativistic corrections. Then for r^^ <C 1 (or M <C Mq) the source 
has size GM/c^ and the flux ^ [A/'i)7i{GM / c^fj^ oc M^''^. For > 1 the source size 
is set by the photosphere radius r^h where aa{r)dr = 1 (i.e. for M ^ Mq, but not so 
large that Qd'^ph) < 0.5 so that our emissivity approximation fails). Then Tph oc M^/^/z/^, 
and the flux oc rl^B^{rph) ^ M^'^li^. 

These simple scaling laws, unfortunately, are not a good description of the variation of 
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source size and flux with M. There are at least three reasons for this. First, relativist ic 
effects are very important; for M ^ Mq emission comes from close to the photon orbit and 
the source size is determined by Doppler beaming and gravitational lensing. Second, for the 
best-bet model at 1.3mm r^^ ^ 1, so in a turbulent flow there is a complicated variation of 
the size of the effective photosphere with M. Third, the emissivity is not precisely frequency 
independent near peak. We therefore need to turn to numerical models. 

The best-best model is taken from a survey of 2D models. Here we adopt the best-bet 
model parameters (a^ = 0.94, i = 85 deR a nd T^/Tg = 3) and use them to s et parameters 
for a fully 3D model ( Shiokawa et al. 2011 ) 0. We use the same data set as Dolence et al 



( 120 12I ). and choose three representative snapshots taken at times when the flow is quiescent 
(t = 5000,9000 and 13000GM/c^ where GM/c^ = 20s). We then recalculate disk images 
and spectra for M = (0.5, 1, 2, 4, 8, 16, 32, 64)Mo. 

The 3D model with M = Mq is broadly consistent with observational data but is slightly 
more luminous at higher energies than 2D models (/3 is lower in the 3D models, and this 
changes the X-ray to millimeter color). The model is self-consistent only for M < 64Mo. 
At higher M the efficiency of the flow is > 0.1 and therefore our neglect of cooling in the 
underlying GRMHD model is not justified. At higher M the 1.3mm photosphere also lies 
outside the limited range in radius where {dM/dr) = 0, so the fiow is not in a steady state. 

The images and total fluxes emitted by the disk at 230 and 345GHz are calculated 



using a ray tracing scheme (INoble et al.l 120071 ). To estimate the size of the emitting region 



we calculate the eigenvalues of the matrix formed by taking the second angular moments of 
the image on the sky (the principal axis lengths). The major and minor axis eigenvalues, 
o"! and (72 respectively, are related by a = FWHM/2.3 to the FWHM of the axisymmetric 
Gaussian model used to interpret the VLB! observations. We use {a) = {ai + cr2)/2 to 
measure the average radius of the emitting spot. 

Figures [1] and [2] show the variation of 230 and 345 GHz Sgr A* model images with M 
(based on a single snapshot from the 3D GRMHD simulation). Evidently modest increases 
M will make the ring-like feature that is the observational signature of the photon orbit easier 
to detect. For M > SMq, however, the ring (or black hole silhouette) is hidden beneath the 
synchrotron photosphere at 230 GHz. The silhouette survives for higher M at 345 GHz, 
disappearing only at M > I6M0. For low M the silhouette is also difficult to detect because 
the emitting region is too small. 



3D GRMHD model parameter survey is still too computationally expensive. It also has a poor return 
on investment given the electron distribution function uncertainties. 
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Fig. 1. — : Horizon silhouette detectability at 230 GHz for various M. Panels from a) to h) 
show the images of Sgr A* calculated for M = (0.5, 1, 2, 4, 8, 16, 32, 64)Mo, respectively. The 
center of the circle is positioned at the image centroid and its radius, r = (ai + 0-2) /2 is the 
RMS radius of the emitting region. 

Figure [31 shows the accretion flow image size ((cr), circle radii in Figured] and [2j) and flux 
at 230 and 345GHz as a function of M. Different point types correspond to three diflFerent 
snapshots from the 3D GRMHD simulation. The size of the emission region and flux increase 
with M . The 345/230 GHz flux ratio increases with the increasing M; this is caused by the 
shift of the synchrotron peak towards higher energies at higher M . 

The following simple fltting formulas describe how a and Fy depend on M: 



230GHz 



M \0.38 



15.2 X (f-) 

_21.1xlogio(f ) + 13, 



for^ < 2 
for# > 2 

Mo — 



[//as] 



(1) 



(c^) 345GHz — 



X 



19.7 xlogio(f) + 8.3, 



for^ < 2 
for^>2 ^ ^ 

Mo — 
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Fig. 2. — : Same as in Figure [T] but for v = 345 GHz. 



S0J2xlog,„(S + 3.8. tor^>2 l'^' 



(3) 



Mo 



M M.3 



^345GHz 



3-2 X (f) 



^0 Jy 

25.13 xlogio(f)+ 0.54, forf > 2 



(4) 



The above constants are nontrivial to interpret because they encapsulate the complexities of 
the accretion flow structure and relativistic eflFects in the radiation transport. The advantage 
of the above formulas is their simplicity. The fltting functions are shown in Figure [3] as dashed 
and dotted lines. 

Figure [4] shows the relation between two observables, a and F^. The size is a linear 
function of the flux and increases more steeply at 230 GHz than at 345 GHz. We also 
provide two phenomenological scaling laws fltted to the data: 



{(y)2mGHz = l-72(F230GHz/Jy) + 9-3 [/xas] 



(5) 
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{(7)MbGHz = 0.73(F345GHz/Jy) + 9-2 [yUas] 



(6) 



Notice that these fits apply to the best-bet model only. For other a*, T^/Tg, or i the scalings 
will be slightly different, although if the source is optically thick the linear scaling follows 
directly from Tg oc 1/r. 

Notice that Equation [5] predicts the change of Sgr A* size by 9 per cent when 230 GHz 
flux changes from 2 to 2.7 Jy. Taking into account observational and theoretical uncertainties 
this is consistent with the observed variations of th e size of the sou rce which increases by a 
few per cent as the flux increases from 2 to 2.7 Jy ( IFish et al.ll201ll ). 



3. Spectra 



Our model spectra are generated by thermal synchrotron emission in the submillimeter/far- 
IR bump, and by Compton scattering in the X-ra ys. The spectral slope of flaring NIR emis- 
sion, and its high degree of linear polarization ( IDodds-Eden et al.ll201ll ). imply that it is 
syn chrotron fro m a small, nonthermal tail of high energy electrons that is not (but can be; 



see 



Leungll2010l ) included in our best-bet model. 



What are the expected scaling laws? Ag am, jjy ' ^ I 
the synchrotron peak is Lp^ak ^ ^TiVpeakju^.aki^^ / ^^f 
that a^GM/c^ = 1. 



n^B ^ M^/^. The luminosity around 
M^/^, where i^peak ^ M^/^ is such 



The emission rightward of the MIR/NIR is produced by Compton up-scattered syn- 
chrotron radiation. The Thomson depth Tgc = gth^^cGM / ^ M, the X-ray luminosity 
is expected to scale as: vLy{v ^ IQ-^^Hz) ^ Lp^ak^sc ^ M-^^/^, assuming that X-rays are 
produced primarily by singly scattered synchrotron photons. 

What do the numerical models show? Figure [5] shows spectra emitted from the 3D disk 
model as obser ved at i ^ 85 deg. T he SEDs are calculated using a general relativistic Monte 
Carlo scheme dPolence et allboOQh . The NIR luminosity iyL^{iy = lO^^Hz) - M^ ^ which 
is only slightly steeper than the expected dependence for the synchrotron peak. vLy{v ^ 
lO^^Hz) ^ M^-^^ agrees well with the expected scaling. 

We conclude that Sgr A* would become a persistent MIR and X-ray source (above the 
present upper limits of 84mJy in MIR and 2.4 x lO^^ergs"^ in X-rays) if M > 2Mo. This is 
conservative in the sense that our models are strictly thermal. The addition of a high energy 
nonthermal tail would only increase the MIR/NIR and X-ray flux. 

We do not consider higher accretion rate models because for M = 64Mo {vLy{v = 
lO^^Hz) = lO^^ergs"^) the model becomes radiatively eflicient, e = Lboi/Mc^ > 0.1, and our 
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neglect of radiative cooling in the underlying 3D GRMHD model is unjustified. 

Finally, notice that the MeV flux increases sharply with M. This suggests that electron- 
positron pair production by photon-photon collisions in the funnel over the pole s of the hole 
would increase sharply (the pair production rigi ^ L^, iMoscibrodzka et al.ll201ll ). If this pair 
production is connected to jet production, it is reasonable to think that a high M Sgr A* 
might also produce a jet. 



4. Discussion 

In summary, we have used general relativistic disk models and relativistic radiative 
transfer to recalculate millimeter images and spectra for Sgr A* at a range of M > Mq. 
Our models predict the following: (1) if the 230 GHz flux increases by more than a factor 
of 2, corresponding to an increase of M by more than a factor of 2, the central accretion 
flow will become a persistent, detectable, MIR and X-ray source; (2) the photon orbit, 
which produces the narrow ring of emission visible in Figures 1 and 2, becomes easier to 
detect for modest increases in the 230 GHz flux; (3) the photon orbit is cloaked beneath 
the synchrotron photosphere at 230 GHz for M > 8Mo, or 230 GHz flux > 13 Jy; (4) the 
photon orbit is cloaked at 345 GHz only at higher M > IGMq, or 230 GHz flux > 17 Jy; 
(5) the size of the source increases in proportion to the flux at both 230 and 345 GHz. We 
suspect that almost any accretion model for Sgr A* with a spatially uniform model for the 
plasma distribution function will reach qualitatively similar conclusions, but that jet models 
may differ signiflcantly. There are order-unity uncertainties in our model predictions due to 
uncertainties in the plasma model. 

What range of changes in M are reasonable? In our best-bet model the mass at radii 
within a factor of two of the pericenter radius is ^ 10~^-^ Mq, assuming steady mass inflow 
from Tp to the event horizon. The addition of even a fraction of the inferred cloud mass 
to the accretion flow in a ring near would (eventually) increase M by a factor of ^ 300. 
On the other hand, stellar winds s upply mas s in t he ne ighborhood of the centra l black 



hole at - 10-^ MQyr-\ Models by lOuataertl (12004h and IShcherbakov fc Baganoflj (12010h 



suggest that most of this mass is ejected in the form of a wind, and that ^ lO~^-^M0?/r~^ 
to 10~^'^MQyr~^ flows inward. A reasonable extrapolation of these models suggest the 
accretion flow at r < has a mass of ^ 2M0; this is comparable to estimates of the mass 
of the inflowing cloud, so in this case we might expect a factor of 2 increase in M. 



This work was supported by the National Science Foundation under grant AST 07-09246 
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Fig. 3. — : Sizes (in terms of a) of the emitting region and the total flux as a function of 
M at 230 (upper panels) and 345GHz (lower panels). Different point types correspond to 
three different snapshots from the 3D GRMHD simulation. The dashed/dotted lines show 
the best fit to the data by Equations [U El |3l and HI 
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Fig. 4. — : Relation between two observables: the flux and the size of the image at v =230 
and 345 GHz. The dashed and dotted hnes are the best flts to the data points. 
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Fig. 5. — : Spectrum emitted by 3D disk model for various M. The Ml Mn is shown on the 



ri^hthand side . Observational points and upper hmits are taken from: iFalcke et al. 



1998. 



An et al.ll2005l . iMarrone et al.ll2006l . iMeha fc Falckell200ll . ISchodel et al.ll201ll.lBaganoff et al 
2003 . The black symbols in the NIR sh owing the flaring stat e are from ICenzel et al" 2003 
An example of X-ray flare is taken from Isaganoff et al.lboOll . 



